Quench dynamics as a probe of quantum criticality 
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Quantum critical points of many-body systems can be characterized by studying response of the 
ground-state wave function to the change of the external parameter, encoded in the ground-state 
fidelity susceptibility. This quantity characterizes the quench dynamics induced by sudden change of 
the parameter. In this framework, I analyze scaling relations concerning the probability of excitation 
and the excitation energy, with the quench amplitude of this parameter. These results are illustrated 
in the case of one-dimensional sine-Gordon model. 



Manipulation of quantum many-body systems, a long 
standing goal of quantum physics, requires understand- 
ing dynamical processes induced by changing external 
parameters of the system. The dynamics is particularly 
interesting and complex when the system is close to crit- 
icality. Recent experiments on cold gases concentrated 
on the many-body dynamics at the Mott-insulator- 
superfluid transition for lattice bosons [TJ, the forma- 
tion of ferromagnetic order in spinor condensates [2], 
and the dynamics of Bose condensation involving sponta- 
neous formation of vortex pairs in the process of thermal 
quench [3 . These experiments illustrate that there are 
several time scales defining the dynamics. 

The dynamics induced by parametric perturbation is 
dictated by the dependence of the Hamiltonian on param- 
eters and also by the structure of excitations for a specific 
system. The relation between the parametric time and 
the intrinsic time scales, characterizing the collective dy- 
namics and the relaxation processes, define the two limit- 
ing regimes of adiabatic (slow) and sudden perturbation. 

In the adiabatic limit, where the parametric time is 
large compared to the intrinsic time scale, the dynam- 
ics of a generic many-body system can be analyzed us- 
ing adiabatic perturbation theory. As a result, one finds 
that upon crossing a quantum critical point the dynam- 
ics is characterized by the rate of change of the external 
parameter and by the critical exponents [4j [5]. Specifi- 
cally, the number of excitations at large times becomes a 
universal power- law functions of the rate. This analysis 
has been recently extended to the non-linear time depen- 
dence of the control parameter [HJ [7], and also to the 
dynamics of open quantum systems [8] . An example of a 
realistic system where this scaling applies is provided by 
the problem of adiabatic loading of strongly-interacting 
one-dimensional bosons into a commensurate optical lat- 
tice [5] that can be mapped to the sine-Gordon model. 

The universality of adiabatic dynamics can be compro- 
mised in realistic situations. First, the relaxation time 
can be comparable or smaller than the parametric time, 
especially in higher-dimensional and non-integrablc sys- 
tems, where thermalization dynamics and also the effects 
of environment [8 complicate the analysis. The details of 
the parametric path at a QCP may also become impor- 
tant, and the scaling results obtained for a single control 



parameter does not apply. 

These limitations are not present for sudden quenches 
of parameters, i.e. fast parametric changes. In this case, 
only the overlap of initial and final states of the system 
are important for characterizing the dynamics. The the- 
oretical analysis of such dynamics has been carried out 
for specific solvable models. Specifically, the ordering dy- 
namics in the systems with spontaneously broken sym- 
metry has been considered in the problem of developing 
superfluidity after crossing the Mott insulator-superfiuid 
transition in the Bose-Hubbard model (TD} E], within 
the BCS model of superconductivity where the dynam- 
ics is induced by abrupt change of the coupling con- 
stant [TH H3] , and also in the problem of magnetization 
ordering in spinor condensates subject to time-dependent 
magnetic field [14] . These examples illustrate sensitivity 
of the dynamics to the initial state and final Hamiltonian. 
Later, the analysis of correlation functions has been ex- 
tended to include spatial dependence in one-dimensional 
integrable models [TBI [S3 CZl HHJ US] which confirmed the 
physical picture of light-cone correlations [TS] . 

In this work, I study the universal scaling in case of 
sudden change of a control parameter, by which the sys- 
tem is driven away from the QCP. Specifically, I consider 
a e?-dimensional quantum system described by a Hamil- 
tonian 7i(A), where A is an external parameter such that 
A = corresponds to the QCP of the system, and the op- 
erator V = d\H is a sum of local operators. For A ^ the 
the system is characterized by finite correlation length, 
£corr ~ \M~ U where v is the critical exponent. The dy- 
namical exponent z relates the critical relaxation time 
to the value of the control parameter, ~ |A| _lyz (see 
Ref. [IB EU 122 for details). 

The main result of this work is the scaling relations 
for the probability of excitation w exc after the control 
parameter is suddenly changed from the critical point 
A c = to a finite value A ^ away from the critical point, 
which is related to the scaling of fidelity susceptibility \f 

w exc « A 2 x F (A)~L rf |Ar«l, 
XF ~ L d \\\ d »- 2 , (1) 

where L is the system size. These results are applicable 
in the vicinity of the critical point, w exc 1, where the 
physical properties are fully characterized by the crit- 
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ical exponents of the transition, and one can rely on 
the perturbative analysis. This scaling is closely related 
to the recently introduced fidelity susceptibility of the 
ground-state \f, that measures sensitivity of the ground- 
state wave-function to the presence of a quantum critical 
point [251123]. 

The main result (JT|) is illustrated in the case of the 
sine-Gordon model at the QCP separating the gapless 
and the gapped phases as a function of the amplitude of 
the symmetry-breaking perturbation (see below for de- 
tails). In this model, I find that fidelity susceptibility 
of the ground-state demonstrates rather peculiar non- 
analyticity: it diverges or vanishes at the QCP depending 
on the coupling constant of the model. This prediction is 
in drastic contrast with the known results regarding the 
behavior of this quantity at the QCP. 

The probability of excitation after a quench can be 
obtained using the perturbation theory. Assuming that 
the final value of the control parameter A ^ is close to 
the critical value A c — 0, we express the initial state in 
the basis of the final states and sum the probabilities of 
the excitation in this basis: 



UK, 



A- 



E 



<0,A|V|n,A) 



E (X) ~ E n (X) 



A 2 Xf(A), (2) 



where we used the expression for fidelity susceptibil- 
ity |24j in terms of the matrix elements of the pertur- 
bation operator. Here, |0, A) and |n, A) are the ground 
state and the excited state, -Eo(A) is the ground state en- 
ergy, E n (X) is the energy of the excited state at a given 
value of A. 

One obtains the scaling of w exc for a sudden quench 
employing the scaling analysis valid close to QCP's [20] 
[5TJ[55]. In this approach, by using the analogy between 
the quantum theory at criticality and the corresponding 
classical field theory in a higher dimension, one obtains 
the scaling dimensions of various quantities. The key ob- 
servation connecting the scaling analysis to the quench 
dynamics is that the scaling transformation is in fact 
equivalent to the dynamical problem of quenching the 
control parameter from the critical point to a finite value 
away from the critical point. According to Ref. [23], the 
scaling relation for fidelity susceptibility away from the 
critical point for the perturbation expressed as a sum of 
the local operators, V = J2 X v ( x ) takes on the form 

J_d\^ |(2Av— 2z— d)v 



XF 



(3) 



where z is the dynamical exponent, and Ay is the scaling 
dimension of V . 

It turns out that the scaling dimension of the pertur- 
bation operator V can be excluded from this relation. 
Indeed, we first notice that the scaling dimension of the 
operator V is obtained from the Hcllmann-Feynman the- 
orem [55] for the derivative of the ground-state energy 



where £qs = Ecs/L d is the density of the ground- 
state energy. Assuming that the scaling dimension of the 
ground state energy is dominated by the perturbation op- 
erator, one obtains Eqs ~ |A| I/Av+1 ~ |A| 2_Q , where a is 
the critical exponent of the ground-state energy [5T] [55] . 
The well-known hyperscaling relation [5TI 155] between 
the exponents, 2 — a = [d + z)u, leads to the identity 



A v = d+ z - l/is. 



(5) 



Substituting this relation into Eq. ([3|, we arrive at the 
main result in Eqs. 0. We notice that the expression for 
the probability of excitation should be much smaller than 
unity, which puts a strict limitation on the final value of 
the control parameter after a quench, that depends on 
the system size. 

The scaling of the excitation energy E exc = 
J2 n [E„(X) - E (X)} w n directly follows from Eq. (j2j, 



where w e 



E 



One obtains 



E„ 



L d \\\ 



(d+z)v 



(6) 



which is identical to the scaling of the ground-state en- 
ergy E GS . 

Similarly, fidelity susceptibility at the quantum criti- 
cal point Xf = Xf(A = 0), scales with the system size 
according to 



Xf 



L 2 '\ 



(7) 



These results are illustrated in the case of the quantum 
sine-Gordon model [55] defined by the Hamiltonian (% — 

1): 



n 



dx[U 2 (x) + (d x (p) 2 -X cos Pip], (8) 



dx£. 



XCGS 



(Q\\V\0\)/L d ~i 



corr 



(4) 



where integration extends over the system size L (the 
boundary conditions are not important for the purposes 
of this work) , and the inter-particle distance defining the 
large-momentum cut-off of the theory is set to unity. The 
canonically conjugate variables II and ip obey the stan- 
dard commutation relation [ip(x), Tl(x')] — i5(x — x'). 
The coupling parameter (3 defines two possible physical 
regimes of the model: the repulsive one (47r < 1 < 8ir) 
in which the excitation spectrum of SG model consists 
of solitons and antisolitons, and the attractive one (0 < 
(3 2 < Att) in which the bound states of soliton-antisoliton 
pairs (breathers) also form. Individual solitons (antisoli- 
tons) have topological charge +1 (—1), while breathers 
have zero topological charge. The spectrum of quasi- 
particles is massive at < f3 2 < 8tt and massless at 
f3 2 > 8tt. Sometimes, it is more convenient to introduce 
the Luttinger parameter K — P 2 /(Att), so that the mas- 
sive regime of the SG model corresponds to < K < 2. 
This model describes e.g. the system of interacting lat- 
tice bosons of the particle density commensurate with 
the lattice spacing [57] , 
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In the following, I analyze the massive regime and 
study the critical properties of the model as one changes 
the control parameter A close to the critical point A c = 0. 
This QCP separates the massive regime of the model at 
A ^ from the massless one at A c = that is the Lut- 
tinger liquid [27]. The dynamical exponent z = 1 as it 
follows from Eq. ([8]), while the critical index v may be 
extracted from the expression for the correlation length 



|A|" 



v= 1/(2- K). 



(9) 



The spectral gap is identified with the soliton mass 
M s ~ £"J. r - lAI 1 /' 2 -^) [H that enters the excitation 
energy of quasi-particles (solitons, antisolitons, and also 
breathers at K < 1) E p = ^Jp 2 + M$. 

The ground-state fidelity susceptibility of the sine- 
Gordon model Q is derived using Eq. ([2|. For the sake 
of clarity, we first discuss this calculation in the repul- 
sive regime (1 < K < 2) of the model. Since the per- 
turbation operator V = — J dx cos (3 tp(x) conserves the 
topological charge and the total momentum of the quasi- 
particles, the dominant contribution to the sum over the 
excited states in Eq. Q is provided by the matrix ele- 
ment (0X\V\p,—p; A) between the ground state |0A) and 
the excited state that contains a soliton-antisoliton pair 
\nX) = \p, —p; A) with zero total momentum 



XF 



(0A|cos f3<p\p, —p; A) 



2E P (X) 



(10) 



where the summation is over half the relative momen- 
tum p of the quasi-particles, and E p = \/p 2 + M 2 is 
the excitation energy of soliton (antisoliton). Since we 
are interested in the scaling of the fidelity susceptibility 
with the control parameter, it is natural to employ the 
well-known result |26] for the scaling dimension (K) of 
operator cos /3ip: 

<OA|COS/^|OA) ~ Corr ~ IM^^K (H) 

where we substituted the correlation length |9]) to obtain 
the scaling with the control parameter. Using this rela- 
tion, it is straightforward to obtain the scaling of the ma- 
trix element involving the soliton-antisoliton pair, which 
at small momentum p <gC £cJrr assumes the form 

(0A|cos/?v?|0A) |A|^/( 2 - A ') 



ALL 



(0A|cos/3<p|p, -p\ A) 

(12) 

Here, the factor M S L takes into account the normaliza- 
tion of the excited state \p, —p; A) , while another factor 
comes from the vacuum average. The large-momentum 
behavior of the matrix element |29j ensures the conver- 



gence of the sum in Eq. ( 10 I for all values of the Luttinger 



parameter < K < 2 in the massive regime of the model. 



Substituting these expressions in Eq. (10 1, one obtains 



XF ~LMj-^_ i|A |(2*-3)/(2-*)_ (13) 



In this relations the first factor, LM S , comes from the 
summation over the states, while another factor comes 



from the matrix element (12) divided by the excitation 
energy. 

It is straightforward to include 2n soliton-antisoliton 
pairs into this calculation (the number is even due to 
conservation of the topological charge, n — 2,3, ..). Due 
to conservation of the total momentum, which is zero in 
the ground state, one needs to sum over only 2n — 1 mo- 
mentum states in Eq. ([2| , since one of the 2n momenta 
is related to the others by the equation Pi + + Pm — 
0. In complete analogy to Eq. (12 1, we obtain the 
scaling of the matrix element | (0| cos/?</?(0)|pi..p2n)| ~ 
\\\ K /( 2 ~ K > /{M s L) n . Upon substitution into the expres- 
sion for the fideleity, one immediately finds the same scal- 



ing (13 1 for all 2n-pair states, as for the one-pair contri- 



bution. We certainly expect, that these multi-pair states 
provide just a numerical correction to the dimensionless 
pref actor of the susceptibility obtained using the one-pair 
state [29]. 

One can extend this calculation to the attractive 
regime at < K < 1, where soliton-antisoliton pairs 
can form bound states, the so-called breathers. The es- 
sential components of the derivation are not affected by 
the presence of the bound states. Specifically, the spec- 
tral gap in this case is also proportional to the soliton 
mass M s ~ lAI 1 ^ 2- ^, and the scaling dimension of the 
perturbation operator in Eq. (11) is also K. Thus, we 



conclude that the scaling of the fidelity susceptibility at 
< K < 1 is also given by Eq. ( 13 1. 



The scaling relation ( 13 1 is in agreement with the gen- 
eral result (JT|), since d — 1 and the critical exponent 
v = 1/(2 — K). The identity ([5| is also satisfied, as it 
follows from the substitution of d = z = 1, the scaling 
dimension of the perturbation operator Ay = K , and 
the critical exponent v. 



Interestingly, upon inspection of Eq. ( 13 1 , we find that 



the character of non-analyticity of the reduced suscepti- 
bility Xf / L at the quantum critical point A c = depends 
on the value of the Luttinger parameter: \f / L diverges 
at < K < 3/2, and it vanishes at 3/2 < K < 2 with A. 

The divergence and vanishing of the fidelity suscepti- 
bility in the corresponding intervals is accompanied by 
the similar behavior of the susceptibility as a function of 
the system size L at the critical point A c = 0, where the 
sine-Gordon model reduces to the Luttinger liquid. In 
this calculation, we first approximate the fidelity suscep- 
tibility by taking into account only the two-phonon con- 



tribution (see Eq. (10 1 where A — 0), and then we study 
the contribution of multi-phonon states. For the Lut- 
tinger liquid, the vacuum average of the operator cos j3(p, 
that defines the behavior of the matrix elements, van- 
ishes in the thermodynamic limit L ^> 1 as dictated by 
the corresponding scaling dimension (K) 



|cos/fy|0) - L 



-K 



(14) 
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The matrix element between the ground state and the 
two-phonon at I > 1 and at small momenta p <!C 1 is 
given by: 



j cos (3<f\p, —p) 



-\p\ 



L\P\ 



\cosf3f\0) 



-\p\ 



-K 



L\P\ 



(15) 

where the factor in the denominator takes into account 
the normalization of the momentum states. The matrix 
element decays exponentially at large momenta p 1 
set by the inter-particle distance (unity in our units). 
Substituting this result in Eq. (JToJ) , where we also use 
E p = p for the excitation energy of the phonons, we find 
that the sum over the states diverges at small momenta 
p c ~ L _1 defined by the system size. The power-counting 
leads to the scaling 



Xf/L ~ L 



3-2K 



(16) 



In analogy with the calculation in the massive phase 
of the sine-Gordon model (A ^ 0), one can include the 
higher-order states in the calculation. Here, the matrix 
element between the multi-phonon and the ground state 
is given by 



Ko|cos#p(o)bi"j>»)l 



-(|pi|+..+|p„|)/2 



-K 



(L\ Pl \..L\ Pn \) 



1/2- 



(17) 



It is straightforward to verify that the dominant contri- 
bution ~ 1? to the fidelity susceptibility is provided by 
the total energy of excitation that enters the denomina- 
tor of Eq. ^ , which after being combined with the other 
factors leads to Eq. (16 1. This result is in agreement 
with Eq. ([7]). We conclude that the reduced suscepti- 
bility chip/L diverges at < K < 3/2 and vanishes at 
3/2 < K < 2 in thermodynamic limit L ^> 1. 

The analysis of the sine-Gordon summarized in 
Eqs. (13 1 and (16 1 demonstrates that the fidelity sus- 



ceptibility is not necessarily a singular function of the 
control parameter at a QCP, although it is certainly a 
non-analytic one. The divergence may be present for 
some values of the critical exponent and absent for other 
ones. Specifically, I find that in the sine-Gordon model, 
the susceptibility diverges at a QCP separating the gap- 
less phase (Luttinger liquid) from the gapped (massive) 
phase, as a function of the amplitude of relevant op- 
erator cos f3 if, when < K < 3/2 and it vanishes at 
3/2 < K < 2. 

In conclusion, I derived scaling relations for the proba- 
bility of excitation and the excitation energy in a system 
driven away from a QCP, and explored its relation to fi- 
delity susceptibility of the ground state, as summarized 
111 Eq. 0. These general results are illustrated in the 
case of sine-Gordon model, in which I find that fidelity 
susceptibility diverges a,t < K < 3/2 and vanishes at 
3/2 < K < 2 with the control parameter, demonstrating 
peculiar type of non-analyticity at the quantum critical 



point. These results can be used to probe the proper- 
ties of QCP's in many-body systems actively studied in 
current experiments on cold atoms. 

I acknowledge discussions with C. De Grandi, V. Grit- 
sev, and A. Polkovnikov at the early stages of this work. 
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